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Abstract 



The ALICE Collaboration has studied J/\// production in pp collisions at ^ = 7 TeV at the LHC 
through its muon pair decay. The polar and azimuthal angle distributions of the decay muons were 
'nI" ■ measured, and results on the I/i/a polarization parameters Xe and Xij, were obtained. The study 

was performed in the kinematic region 2.5 < y < 4, 2 < pi < ^ GeV/c, in the helicity and Collins- 
Soper reference frames. In both frames, the polarization parameters are compatible with zero, within 
^ ■ uncertainties. 
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Almost forty years after its discovery, heavy quarkonium still represents a challenging testing ground 
for models |1] based on Quantum Chromodynamics (QCD). Results obtained for charmonium produc- 
tion at the Tevatron collider in the nineties IS led theory to recognize the role of intermediate quark- 
antiquark color octet states in the production process, in the framework of the Non-Relativistic QCD 
model (NRQCD) [3|. This approach brought the calculations of pt spectra to agree rather well with the 
data [4 1 {pt is the transvere momentum, i.e. the momentum component perpendicular to the colliding 
beam direction). However, the same calculations were not able to reproduce satisfactorily the polariza- 
tion results for the J/y obtained by the CDF experiment at ^/s = 1.96 TeV [5]. In particular NRQCD 
at leading order (LO) predicts for high-/7t 3/W (Pt ^ ^ significant transverse polarization, i.e. a 

dominant angular momentum component 7^. = ±1, the z-axis being defined by the J/y/ own momentum 
direction in the center of mass frame of the pp (pp) collision. Contrary to this expectation, CDF data ll5l 
rather exhibit a mild longitudinal polarization {J^ = 0). In a recent renaissance of quarkonium studies, 
also related to the publication of results from RHIC at y/s = 0.2 TeV [6], next-to-leading order (NLO) 
corrections for both color singlet and color octet intermediate states were calculated, and their impact 
on the Pi spectra was found to be quite important dllllll. The influence of these corrections on the 
polarization calculations is expected to be significant ifTOl ITTl and still has not been completely worked 
out. The start-up of the LHC provides the possibility to perform charmonium measurements in a new 
energy domain, over large ranges in p^ and rapidity {y = 0.51n[(£' + Pz)/{E — Pz)], where E is the en- 
ergy and pz is the momentum component parallel to the colliding beam direction). Various theoretical 
approaches |El[T2l[T3| proved to be rather successful in describing the first LHC experimental results on 
the J/i//^ Pi spectra fT4l lT5l[T6llT7 1. The measurement of polarization clearly represents a more stringent 
test of the theoretical calculations, offering therefore the possibility of confirming/ruling out the current 
QCD approach to charmonium production. 

In this Letter we present the results of a study of J/yA polarization at the LHC, carried out by the ALICE 
experiment in pp collisions at y/s = 7 TeV. The ALICE experiment ifTSll is based on a central barrel, 
covering the pseudorapidity region \^\< 0.9 [19 1 and a muon spectrometer, with 2.5 < T] < 4 coverage. 
The polarization results presented in this Letter refer to inclusive J/i/A, measured via the J/i/a — )• pL^pi^ 
decay in the muon spectrometer. The spectrometer ifTTl consists of a 10 interaction length (Ai) thick 
front absorber, to remove hadrons, followed by a 3 T-m dipole magnet. Charged particles which exit 
the front absorber are tracked in a detector system made up of five stations, each one with two planes of 
Cathode Pad Chambers. The tracking system is followed by a 7.2 X\ iron wall, which absorbs secondary 
hadrons escaping the front absorber and low-momentum muons. Finally, a trigger system, based on 
Resistive Plate Chambers, is used to select candidate muons with a transverse momentum larger than a 
given programmable threshold. 

The analysis presented in this Letter was carried out on a significant fraction of the 2010 sample of muon- 
triggered events, corresponding to an integrated luminosity Lim ~ 100 nb^^. The usual event selection 
cuts, already applied to previous analysis of J/i/a production [17], were also used for the polarization 
study. Events with at least one vertex reconstructed in the Inner Tracker System (ITS) |[20ll are retained 
for the following analysis if they contain at least two tracks reconstructed in the muon spectrometer, out 
of which at least one has to satisfy the trigger condition (1 GeV/c pi threshold). We note that with this 
requirement the acceptance of the spectrometer for J/y/ extends down to pi = 0. The tracks must satisfy 
the condition 2.5 < Tj < 4 and must also have 17.6 < /?abs < 88.9 cm, where /?abs is the radial distance 
of the track from the beam axis at the exit of the front absorber (z = 503 cm). The latter requirement 
eliminates forward tracks which, due to the high-Z material used in the absorber in that region, are 
strongly affected by multiple scattering. Finally, a rapidity cut 2.5 < j < 4 is applied to the selected 
muon pairs. 
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The distribution of the J/i/a decay products can be expressed in its general form [21] as 

oc \ (1+Aecos^e + A0sin^ecos20 
3 + Ae 

+Ae0 sin20 cos (j)) , (1) 

where 6 (0) are the polar (azimuthal) angles in a given reference frame. In this analysis, the CoUins- 
Soper (CS) and helicity (HE) frames were considered. In the CS frame the z-axis is defined as the bisector 
of the angle between the direction of one beam and the opposite of the direction of the other one, in the 
rest frame of the decaying particle. In the HE reference frame the z-axis is given by the direction of the 
decaying particle in the center of mass frame of the collision. The = plane is the one containing 
the two beams, in the J/y rest frame. Equation [T] contains the three parameters Xg, and Xq^, which 
quantify the degree of polarization. In particular, Xq > Q values indicate transverse polarization, while 
a longitudinal polarization gives Xq < 0. In principle the values of the parameters could be extracted 
by means of a fit to the acceptance-corrected 2-D distributions for cos 6 vs (j). However, the limited 
statistics (about 6.8 • 10^ signal events in the pt range under study) makes a 2-D binning impossible. 
Therefore the study of the angular distributions was separately performed on the polar and azimuthal 
variables. In particular, Xq and were obtained by studying the distributions 

W(cose)oc (i + Aecos^e) 

3 + Xe 

2X 

1 + — ^cos2(^», (2) 
obtained by integrating Eq.[T]in the ^ and cos 6 variables, respectively. 

The distributions of the angular variables for the J/y decay products were obtained starting from the 
study of the dimuon invariant mass spectra. The study was performed in five bins for the | cos 6 \ variable 
(the angular distribution is symmetric with respect to cos 6 = 0), in the range < | cos 0| < 0.8. For the 
azimuthal variable four bins in |(^| were defined, in the range < \ ^\ < n/2 (values between 7l/2 and n 
were mirrored around \ = n/2, due to the period of the cos20 function). The analysis was carried out 
in three transverse momentum intervals (2 < pi <3 GeV/c, 3 < pt < 4 GeV/c and 4 < pt < S GeV/c). 
The limits of the explored pi range are related to the strong decrease of the acceptance for large | cos 6 \ 
values at low pt and to the limited statistics at high pi. 

The number of J/i/A signal events for the various bins in | cos 6 \ and | <p \ were obtained by means of fits 
to the corresponding dimuon invariant mass spectra performed in the range 1.5 < m^^ < 5 GeV/c^, and 
in Fig. [T] we show one of them as an example. The J/i/A signal was described by a Crystal Ball function 
(CB) f22| while for the background an empirical function, corresponding to a Gaussian with a width 
linearly depending on mass, was adopted. The position of the CB peak was left as a free parameter in the 
fits, and was found to correspond to the nominal J/y pole mass within at most 1%. The width of the CB 
function obtained from the data (between 72 and 120 MeV/c^, depending on the kinematics) was found 
to be in agreement with the Monte Carlo (MC) within ~ 8 - 10 MeV/c^. In the fits, the width of the CB 
function for each bin / (where / represents a certain | cos 6 \ or \ interval for the J/y Pt bin under study) 
was fixed to a^, = aj/,^ • {(^^^ /(^fj^)^ i-C- by scaling the measured width for the angle-integrated 
spectrum with the MC ratio between the widths for the bin / and for the integrated spectrum. The quality 
of all the fits is satisfactory, with ;^^/ndf in a range between 0.63 and 1.34. Signal over background ratios 
in a ib3a mass window around the CB peak vary between 0.5 and 3.5. The number of signal events per 
bin ranges from ~100 (for 2<pt<3 GeV/c, 0.6 < | cos 0cs| < 0.8) to ~1000 (for 2<pt<3 GeV/c, 
0< |cos0cs| <0.15). 

The polarization parameters for the J/\j/ were obtained by correcting the number of signal events Nj^^ for 
each bin for the product A, x £, of acceptance times detection efficiency, calculated via MC simulation. 
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Fig. 1: The dimuon invariant mass spectrum for 2 < pt < 3 GeV/c, < | cos Que \ < 0. 15, together with the result 
of the fit. The conttibutions of the signal and background are also shown as dashed lines. 

and then fitting the corrected angular distributions with the functions shown in Eq. |2] The simulation 
includes, for the tracking chambers, a map of dead channels and the residual misalignment of the de- 
tection elements and, for the trigger chambers, an evaluation of their efficiency based on data. It also 
includes a random misalignment of the tracking detector elements, of the same size of the resolution 
obtained by the offline alignment procedure ifTTl . For both tracking and triggering detectors, the time 
variation of the efficiencies during the data taking period was accounted for (see [17] for details). Since 
the cos 6- and -acceptances are strongly correlated, the acceptance values as a function of one variable 
strongly depend on the input distribution used for the other variable. Given the fact that the correct input 
distributions are not known a priori, but rather represent the outcome of the data analysis, an iterative 
procedure was followed in order to determine them. In the first iteration a flat distribution of the angular 
variables (equivalent to a totally unpolarized J/y/ distribution) was adopted to calculate the acceptances. 
After correcting the signal with those acceptances, a first determination of the polarization parameters 
is performed, and the results are then used in a second determination of the acceptance values. The 
procedure is then repeated until convergence is reached, i.e. the extracted polarization parameters do 
not vary by more than 0.005 between two successive iterations. This occurs, for this analysis, after at 
most three steps. It was also checked that using polarized MC input distributions in the first iteration the 
procedure converges towards the same results as in the default, unpolarized, case. Typical A, x £, values 
vary between ~0.22 (0.05) at low pt and large | cos 0| and ~0.41 (0.63) at large pt and small | cos 0| for 
the HE (CS) frame. 

A simultaneous study of the J/y polarization variables in several reference frames, as first carried out 
in hadroproduction studies by the HERA-B experiment [23], is particularly interesting since consistency 
checks on the results can be performed, using combinations of the polarization parameters which are 
frame-invariant. In particular we made use of the invariant F = {Xg + 3X,j,)/{l — A^) [21 ], performing a 
simultaneous fit of the | cos 6 \ and | (p \ distributions in the two reference systems and further constraining 
the fit by imposing F to be the same in the CS and HE frames. In Fig. |2]we present, as an example, the 
result of such a fit relative to the last iteration of the A, x e,- calculation, for 2 < /Jt < 3 GeV/c. The x^/ ndf 
values {ndf = 10) are 1.08, 1.00, 1.32 for 2 < /^t < 3, 3 < /?t < 4 and 4 < /^t < 8 GeV/c, respectively, 
showing that the quality of the fits is good. Compatible results are obtained when the constraint on F is 
released. 

In the analysis described so far, the Xq^ pai^ameter was implicitly assumed to be zero in the iterative 
acceptance calculation. In the one-dimensional approach followed in this analysis Ae^ could be estimated 
from the data, defining an ad-hoc variable which is a function of cos Q and and contains Xq^ as a 
parameter (see [21] for details). In principle, the iterative procedure applied to Xq and X^ determination 
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Fig. 2: The acceptance corrected angular distributions for the J/y/ decay muons, for 2 < pt < 3 GeV/c. The 
simultaneous fit to the results in the CS and HE frames is also shown. The plotted errors are purely statistical. The 
horizontal bars represent the bin width. 



could be extended to include Xq^; however, in some cases, relatively small statistical fluctuations in the 
distributions of the measured variables tend to induce large variations of the fitted values in the following 
iterations, leading to convergence problems. A check of the Xg^ =0 assumption was done a posteriori for 
each pt bin, by fitting the distributions, corrected with an acceptance which makes use of the measured 
Ae and values as inputs. In this way we get for all the pt bins Xgtj, values compatible with zero for 
both CS and HE reference frames. We also note that all the previous experiments assumed = in 
their analysis, with the exception of HERA-B [i23il who measured it in pA collisions at ^/s = 41.6 GeV 
and found values ranging from to 0.05. 

Various sources of systematic uncertainty on the measurement of the polarization parameters have been 
investigated. The uncertainty on the signal extraction was studied by leaving in the fits the width of the 
CB function as a free parameter. This choice leads to an absolute variation of the polarization parameters 
between 0.02 and 0.10. Another sizeable source of systematic uncertainty is the choice of the input 
distributions for pt and y in the simulation. It was evaluated by comparing the results obtained with a 
parameterization of our 7 TeV results on differential J/y cross sections lilTl with those obtained using 
an extrapolation of lower energy results ||24|| . The absolute effect on the polarization parameters varies 
between 0.01 and 0.07. For the lowest pt bin, the acceptance in the HE frame drops by about 40% in the 
highest I cos 0| bin used in the analysis (0.6 < | cos 0| < 0.8), and has also a strong variation inside the bin 
itself. We therefore followed an alternative approach, fitting the angular spectrum in the restricted interval 
< I cos 0| < 0.6 (instead of the default choice < | cos 0| < 0.8) and we conservatively considered the 
variation in the result of the fit (0.15) as an additional systematic uncertainty on Xq. For consistency, the 
same evaluation was performed in the CS frame. The role of the systematic uncertainties on the trigger 
and tracking efficiency [17 ] was also studied. The first was evaluated by varying the efficiency values for 
each detector element by 2% with respect to the default values in the simulation. This choice is related to 
the estimated uncertainty on the detector efficiency calculation. For the second we have used the rather 
conservative choice of comparing the reference results, obtained with realistic dead channel maps, with 
those relative to an ideal detector set-up. The result is typically 0.03-0.04. Finally, by quadratically 
combining the results for the various sources, values between 0.04 and 0.21 are obtained for the global 
systematic uncertainties. 



J/Y polarization in pp collisions at ^/s=^ TeV 



7 



In Fig. [3] we show the results on Xg and for inclusive J/i/a production. In both frames all the parameters 
are compatible with zero, with a possible hint for a longitudinal polarization at low pi (at a 1.6a level) 
in the HE frame. The numerical values are given in Table [T] 
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Fig. 3: Xg and X,j) as a function of pt for inclusive J/\//, measured in the HE (closed squares) and CS (open circles) 
frames. The error bars represent statistical errors, while systematic uncertainties are shown as boxes. 



Table 1: The values of Xg and X,j) in the two reference frames. Statistical and systematic uncertainties are quoted 
separately. 

PtiiPt)) , , 

(GeV/c) 



2-3 (2.5) 
HE 3-4 (3.4) 
4-8 (5.1) 



CS 



-0.36 ±0.09 ±0.21 
-0.20±0.11±0.13 
0.00 ±0.10 ±0.10 



0.05 ±0.04 ±0.04 
0.01 ±0.05 ±0.05 
0.00 ±0.04 ±0.04 



2- 3 (2.5) 

3- 4 (3.4) 

4- 8 (5.1) 



-0.10±0.14±0.13 
-0.06 ±0.14 ±0.07 
-0.09 ±0.10 ±0.08 



-0.04 ±0.08 ±0.07 
-0.03 ±0.08 ±0.05 
0.03 ±0.06 ±0.07 



The inclusive J/\j/ yield is composed of a "prompt" component (direct J/y + decay of the V^(2S) and 
Xc resonances) and of a component from B-meson decays. In the pi range accessed in this analysis 
the B-meson decay component accounts for 10% (2 < pj < 3 GeV/c), 12% (3 < pj < 4 GeV/c) and 
15% (4 < Pt < 8 GeV/c) of the inclusive yield, according to the LHCb measurements carried out in 
our same kinematical domain ifTSl . The polarization of the non-prompt component is expected to be 
quite small. In fact, even if a sizeable polarization were observed when the polarization axis refers to the 
B-meson direction [25 J . it would be strongly smeared when it is calculated with respect to the direction 
of the decay J/i/a ifTSl . as observed by CDF, who measured in this way Xe{J /y ^ ~ —0.1 in the 
HE frame Q. Assuming conservatively \Xe{J /y ^ B)\ < 0.2 for both frames, and taking into account 
the fraction of the inclusive yield coming from B-meson decays [T5l, the difference between prompt and 
inclusive J/i/a polarization was estimated and found to be at most 0.05, a value smaller than the systematic 
uncertainties of our measurements. Concerning higher-mass charmonia, the Xc — ^ J/v^± / decay cannot 
be reconstructed in the muon spectrometer, and the '^(25) jxpi statistics is currently too low. Values of 
the feed-down ratios measured mainly by lower energy experiments range from ~10% for the y{^^) Il26l 
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to 25-30% for the Xc [27], implying that there could be a sizeable difference between direct and prompt 
J/i/A polarization. 

The results presented in Fig. |3] extend the study of the J/i/a polarization to LHC energies and therefore 
open up a new testing ground for theoretical models. At present, NLO calculations for direct J/i/a po- 
larization at the LHC via the color-singlet channel I1T21 [TOl predict a large longitudinal polarization in 
the HE frame {Xq ~ —0.6) at pt ~ 5 GeV/c, which is in contrast with the vanishing polarization that 
we observe in such a transverse momentum region. The contribution of the S-wave color-octet channels 
was also worked out [9] and indicates a significantly different trend (large transverse polarization) with 
respect to the color-singlet contribution, but again in contrast with our result. In this situation, a rigorous 
treatment on the theory side of all the color-octet terms (including P-wave contributions) is mandatory, 
as well as a study of the contribution of Xc and ^{2S) feed-down which, as outlined before, is important 
for a quantitative comparison with our result |[28ll . Such studies are presently in progress and the com- 
parison of their outcome with the results presented in this Letter will allow a very significant test of the 
understanding of the heavy-quarkonium production mechanisms in QCD-based models. 

In summary, we have measured the polarization parameters Xq and for inclusive J/i/A production 
in y/s = 7 TeV pp collisions at the LHC. The measurement was carried out in the kinematical region 
2.5 <3'<4, 2<pt<8 GeV/c. The polarization parameters Xq and X^ are consistent with zero, in both 
the helicity and ColUns-Soper reference frames. These results can be used as a stringent constraint on 
the commonly adopted QCD framework for heavy quarkonium production. 
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